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Abstract 

The magneto-optical spectra of ABC-stacked trilayer graphene are enriched by 
the electric field. A lot of prominent absorption peaks, which arise from the inter- 
Lnadau-level transitions, gradually change from the twin-peak structures into the 
double-peak ones with the increasing electric-field strength. This comes from the de¬ 
struction in mirror symmetry of xy-plane and the non-equivalence for two sublattices 
with the identical projections. Specially, a single threshold peak becomes a double¬ 
peak structure, owing to the Fermi-Dirac distribution. The absorption frequencies 
continuously grow or decline except for the anti-crossings of Landau levels. Also, such 
anti-crossings can induce extra double-peak structures. 
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Graphene, with a fascinating hexagonal lattice structure, was first fabricated by me¬ 
chanical exfoliation in 2004. [T] Its unique electronic and optical properties can be easily 
modulated by changing the layer number [2], the stacking configuration [3], the external 
magnetic |3] and electric held [5], and the deformed structure pa]. These geometric varieties 
and tuning effects have attracted numerous theoretical and experimental studies especially 
for the magneto-electronic and -optical properties in few-layer graphenes. [Mj The quan¬ 
tized Landau levels (LLs) are strongly dependent on the stacking configuration. An N-layer 
graphene has N groups of conduction and valence LLs. The three highly symmetric layer 
configurations, AAA-, ABA-, and ABC-stacked trilayer graphene, exhibit the totally dif¬ 
ferent magneto-optical properties. The three categories of excitation channels, belonging 
to intra-group optical transistions, are revealed in the AA-stacked trilayer graphene. The 
inter-group LL excitations are absent because of the special relationships among the LL 
wave functions. [10] The ABA trilayer one is expected to present five categories of absorp¬ 
tion peaks since their LLs can be regarded as the superposition of those of the AB-stacked 
bilayer and monolayer systems. Moreover, the ABC-stacked trilayer graphene possesses 
nine categories of absorption peaks, including all the intra-group and inter-group LL exci¬ 
tations. m A specific magneto-optical selection rule, An = 1 in the variation of quantum 
number, is identified in the above-mentioned systems. In particular, there are three kinds 
of the selection rules, An =0, 1, and 2, in the sliding graphene with various stacking config¬ 
urations. [12] On the experimental side, only the intra-group excitation channels from the 
first LL group in the AB-stacked few-layer graphenes (N=l ~ 5) have been verified until 
now. [13] In this work, we focus on how a uniform perpendicular electric held (E = E 0 z) 
affects the low-frequency magneto-optical properties in the ABC-stacked trilayer graphene, 
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including the structure, intensity and frequency of absorption peaks. 

We have developed the generalized tight-binding model for various external fields m 
ns], in which the Hamiltonian matrix is built from the tight-binding basis functions, i.e. the 
subenvelope functions on the different sublattices. It can be further utilized to investigate 
the magneto-optical properties of ABC-stacked trilayer graphene in E. Detailed investiga¬ 
tions have been performed on how the inter-LL optical transitions are diversified by the 
relationship between the interlayer atomic interactions and the electric field. This study 
shows that the electric field causes the splitting of Landau levels, and hence the abundant 
absorption spectra are induced, such as the E 0 -dependent peak structure, frequency and 
intensity. The twin-peak structure or a single-peak structure will be replaced by a double¬ 
peak structure at a sufficiently large electric-field strength. The extra peaks are generated 
by the intragroup LL anticrossings during the variation of E 0 ; futhermore, nduced by the 
abnormal E 0 -dependent absorption frequencies are obtained. The theoretical predictions 
could be verified by the optical spectroscopy methods. pdH ITTT| 

The magnetic held, electric held and interlayer atomic interactions are simultaneously 
taken into account in the generalized tight-binding model. The primitive unit cell of the 
ABC trilayer graphene consists of six carbon atoms (A 1 , B 1 , A 2 , B 2 , A 3 ; B 3 ), where 
superscripts 1, 2, and 3 correspond to the hrst, second, and third layer. The atomic 
interactions include one intralayer type (/3 0 = 3.16 eV) and hve interlayer types (/3 0 = 3.16 
eV, /3i = 0.36 eV, (3 2 — —0.01 eV, /3 3 = 0.32 eV, /? 4 = 0.03 eV; /3 5 = 0.013 eV). The 
vector potential caused by a uniform perpendicular magnetic held introduces a periodic 
phase in the atomic interactions. Thus, the unit cell is enlarged by 2 Rb times (Rb = 
79000/I?o), and the Hamiltonian is a 12 Rb x 12 Rb Hermitian matrix. Moreover, a uniform 
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perpendicular electric field leads to a potential energy difference (V g = E 0 d ; interlayer 
distance d=3.35 A) between two neighboring layers, i.e., V g changes the site energies of 
the diagonal Hamiltonian matrix elements. The eigenvalues and eigenfunctions can be 
solved effectively by a special diagonalizing method which is used to obtain a band-like 
Hamiltonian matrix. 

The optical absorption function according to the Fermi golden rule at zero temperature 
is directly obtained through 
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where the superscripts c and v represent, respectively, a conduction and a valence state, 
f(E c,v (k,n)) is the Fermi-Dirac distribution function, and T (~ 1 meV) is the broadening 
parameter. The absorption spectrum accounts for the inter-LL transitions from the occu¬ 
pied state to the unoccupied one. Only vertical transitions of Ak x = 0 and Ak y = 0 are 
allowed beacuse of the zero photon momentum. The absorption intensity is dominated by 


the velocity matrix, which is calculated by the gradient approximation. 


For the elec¬ 


tric polarization E || y, the velocity matrix elements are obtained from the multiplication 
of the three matrices corresponding to the initial state, the final state and dH/dk y . The 
spectral intensity is dominated by the nearest-neighbor Hamiltonian matrix elements, since 
the in-plane atomic interaction /3 0 makes the largest contribution to the inter-LL excita¬ 
tions. The optical transistions are available only when the subenvelope function of the A 1 
(B l )sublattice in the initial state has the same mode with that of B l ( A *) sublattice in the 
finial state. 

In the presence of an electric field, each LL is twofold degenerate rather than the or¬ 
dinary fourfold degenerate, mainly owing to the destruction in xy-plane mirror symmetry. 
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The six subenvelope functions are localized at 2/6 (5/6; purple) and 4/6 (1/6; blue) posi¬ 
tions, as shown in Fig. 1 for B$ = 48 T and V g = 0.23 eV. The quantum numbers of the 
first LL group, rij g and n c {^ 1 are characterized by the Bj and Aj sublattices, respectively, 
where the subscripts /3 and a stand for the localized positions 2/6 and 4/6. Even with 
the same quantum number, the conduction and the valence LLs have their six subenvelope 
functions behave differently. For example, the amplitudes of the A)- and Aj (Bj and Bj) 
in n v ig = 4 (n c ig = 4) LL are much smaller than those in n c ig = 4 (n v ig = 4) one. In general, 
the subenvelope functions corresponding to the identical planar projections, (Aj -H- Bj) 
and (Aj -h- Bj), are similar in their node numbers and amplitudes. However, the potential 
difference will destroy the stacking symmetry between these two sublattices, so that the 
subenvelope functions are no longer equivalent for a sufficiently large V g . The sharp con¬ 
trast of the amplitude in subenvelope functions is expected to induce dramatic changes in 
the absorption peaks. 

The electric field and the atomic interactions lead to a feature-rich energy spectrum, 
as shown in Fig. 2(a). Some LL energies ( E c,v ) monotonously increase or decrease with 
the rising V gi and the others demonstrate non-monotonic V g - dependences. In particular, 
the intragroup LL anticrossings occur frequently during the evolution of the V/-dependent 
LL energies (black circles), since some LLs possess multi-mode at certain V g ranges. The 
amplitudes of main and side modes in the subenvelope functions have a drastic change. [T2] 
For example, the subenvelope functions, (n c la = 1 & = 4) and (n\ g = 1 & n\ g = 4), 

are strongly hybridized around the anticrossing point (V/ ~ 0.19 eV), the progressions 
are, respectively, illustrated in Fig. 2(b) and 2(c) for 0.15 eV < V g < 0.25 eV. After 
that, their quantum numbers are exchanged. The modes of the subenvelope wavefunctions 
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become well-defined again in the further increase of V g . These variations in the amplitude, 
waveform, and zero-point number of the subenvelope function, which are modulated by V g , 
directly reflect on the absorption spectra. 

The low-frequency absorption spectra with respect to the first group of LLs exhibit 
many more special peak structures as being diversified by the electric field. The threshold 
inter-LL excitation, coming from the optical transition: the occupied = 3 LL to the 
unoccupied n\ = 2 LL, exhibits a single-peak structure for V g — 0, which is attributed 
to the Fermi-Dirac distribution and the selection rule of An = 1. mi The LLs with the 
different localization positions make the same contribution. The LL degeneracy is destroyed 
by V g , therefore, the threshold excitation becomes a double-peak structure. The lower- and 
the higher-frequency comes from n\o = 2 — > n\p = 3 and n\ a = 3 —> n\ a = 2 (arrows in 
Fig. 2(a)) with their subenvelope functions localized at 2/6 and 4/6 positions, respectively. 
The peak structure stays unchanged, and their frequencies increase with the increment of 
V g . However, the other absorption peaks for V g = 0 possess the twin-peak structures with 
the similar absorption intensity, HH which is caused by the asymmetric energy spectrum 
between the occupied and unoccupied LLs. Such structures also emerge in two different 
localiztion positions, 2/6 and 4/6 for a small V g . As denoted by (3^4/3 & 4^3^) and (4 Q 3 Q & 
3 q 4 q ) in Fig. 3, two twin-peak structures or four absorption peaks are displayed at V g = 0.01 
eV. With the further increasing V g , they are replaced by a double-peak structure (black 
circle at V g = 0.11 eV) since a certain peak of the twin-peak structure is almot vanishing. In 
other words, the two peaks in a double-peak structure, respectively, arise from the inter-LL 
transitions with the distinct localization centers of 2/6 and 4/6. Apparently, the dramatic 
changes in peak structure indicate the strong competition between the interlayer atomic 
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interactions and the electric field. 


The V^-induced rich absorption spectra deserve a closer investigation. For the equivalent 
sublattices, (A] O Bj) and (A? O B : - ), with the identical planar projections possess the 
similar waveforms and amplitudes in the subenvelope functions. m The equivalence will 
be destroyed when a electronic static gate voltage induces a potential energy difference on 
the distinct layers. At a sufficiently large V g , the amplitudes of A 1 and A 2 in the valence 
LLs and those of B 2 and B 3 in the conduction LLs become very small (red rectangles in 
Fig. 1). And the spectral intensity of the absorption peak 3^4^ (4 Q 3 Q ) is proportional 
to the inner product calculated from the B l sublattices of rip = 3 (n v a = 4) and the A 1 
sublattices of njj = 4 {n c a = 3). However, for the absorption peak 4^3^ (3 Q 4 Q ), the intensity 
is determined by the inner product of the A 1 sublattices in n v o = 4 ( n v a = 3) and the B l 
sublattices in njj = 3 ( n c a = 4). As a result, the latter is much weaker than the former. 

The electric field induces the intragroup LL anticrossings and thus extra double-peak 
structures. The existence of extra absorption peaks is mainly controlled by the weight 
of the side mode in the progressive LL. For example, the spectral intensities of lphp and 
5 a l a (brown circles in Fig. 3) begin their gradual rises from V g — 0.15 eV and reach their 
maximums when V g gets to 0.19 eV, since the subenvelope functions of n =1 and 
=4 are strongly hybridized with each other (Fig. 2(b) and (c)). After that, the intensities 
drop low under the increasing V g and disappear when the quantum numbers of subenvelope 
functions restore to be well-defined again. 

The low-lying inter-LL absorption frequencies (conn) strongly depend on the electric-field 
strength. The absorption frequencies of the double-peak structures monotonously decrease 
as V g increases, while the opposite is true for the threshold double-peak ones, as shown in 
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Fig. 4. Such dependences are revealed to originate from the V^-dependent LL energies. 
The perpendicular external electric field induces the opening of energy gap; therefore, the 
absolute energies of the six valence and conduction LLs respectively located at the 2/6 and 
4/6 positions, which are near the Fermi level and designated to the quantum numbers 0, 
1 and 2, will rise with the increasing V g . However, the behavior of the others LLs at the 
2/6 (4/6) position is just opposite to these six LLs. As a result, the energy differences of 
the LLs can be clearly reflected on uj Rn . Moreover, a;„ n exhibits some discontinuous In¬ 
dependent structures at a certain electric field strength (black circles), mainly owing to the 
intragroup LL anticrossing. The aforementioned characteristics of the V^-induced magneto¬ 
absorption spectra, like the peak structure, peak intensity, and V/-dependent absorption 
frequency, could be verified by optical spectroscopy methods, such as transmission [2TLI22] , 
reflection J23j[2l] and Raman scattering spectroscopies. [16^25,126] 

The generalized tight-binding model is developed to investigate optical properties of 
ABC-stacked trilayer graphene in magnetic and electric fields. The prominent LL magneto¬ 
absorption peaks are diversified by the electric field. With the increased V g , the twin-peak 
structures gradually change into the double-peak ones. The former and the latter, which 
correspond to the inter-LL transitions with the same and the different localization centers, 
are, respectively, caused by the asymmetric LL energy spectrum and the destructions in 
both the mirror symmetry and equivalence of two sublattices. The special peak structures 
are mainly determined by the competition between the interlayer atomic interactions and 
the electric-held strength. Especially, a threshold single peak is replaced by a double-peak 
structure because of the Fermi-Dirac distribution. Moreover, the anti-crossings of LLs can 
create extra double-peak structures and thus disrupt the dependence of peak frequency on 


Vg. The experimental examinations of the V^-induced magneto-optical spectra are useful 
in identifying the stacking configurations among few-layered graphene systems. 
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Figure 1. The first group of Landau levels and their respective six subenvelope 
wavefunctions for Bq = 48 T and V g = 0.23 eV. 
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Figure 2. (a) The two Landau-level subgroup spectra at B 0 = 48 T as the 

function of the electric field (V^,) illustrated in two different color codes: purple 
and blue. The a-type and /3-type intragroup Landau-level anticrossings are 
indicated by red and green rectangles, respectively. (b)The variations of the 
Aj subenvelope wavefunctions for n c a — 1 and n c a = 4 with respect to the a-type 
intragroup Landau-level anticrossing. (c)The variations of the Bj subenvelope 
wavefunctions for rip = 1 and n /3 — 4 with respect to the /3-type intragroup 
Landau-level anticrossing. 
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Figure 3. The progression of the absorption peaks in the low-frequency region 
from V g = 0.01 to 0.23 eV. 
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Figure 4. The V^-dependent frequencies of the low-lying inter-Landau level 
absorption peaks. 
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